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Abstract: We analyze the physical meaning of fluctuations of the phase space con- 
traction rate, that we also call entropy creation rate, and its observability in space-time 
intermittency phenomena. For concreteness we consider a Navier-Stokes fluid. 

§1. The chaotic hypothesis in turbulence. 

Consider a Navier-Stokes (NS) fluid in a container V which we take, for simplicity, cubic 
with periodic boundary conditions and subject to a constant volume force f ^ {x) with 
max I V' ( 2; ) I = 1 and with only Fourier harmonics corresponding to large wavelength of 
the order of the linear size L of V . 

The viscosity will be denoted but it is convenient to rescale space, time, velocity and 
pressure to write the equations in dimensionless form as 

u+Ru-du^Au-dp+ip, d-u = 0, R = JL^v^'^ (1.1) 

in a container of side L — 1, where R is the Reynolds number. We can suppose that 
Jy u dx = (because of translation invariance). 
We assume the chaotic hypothesis 

Chaotic hypothesis: Asymptotic motions of a turbulent flow develop on an attracting 
set A in phase space on which time evolution u ^ St u can be regarded as a transitive 
Anosov system for the purposes of computing time averages in stationary states. 

Here we investigate some assumptions under which the hypothesis acquires some non 
trivial predictive value with implications that can have experimental relevance. For 
earlier reviews on the chaotic hypothesis see [Ga98a] , [Ga96d] , [Ga99a] . A recent one is 
[Ru99a]. 

§2. The OK41 cut off. 

Anxiety often mars the beginning of any discussion on the NS equations: it is a fact that 
to date there is no theory that allows a constructive solution of the equations via a con- 
trolled approximation scheme. Nevertheless most people rapidly recover and adopt the 
viewpoint that "physically there is an effective ultraviolet cut-off" and the NS equations 
can be reduced to ordinary equations: 

The OK41 cut— off hypothesis: There exists kq > such that if the NS equation, 
(1-1), is truncated in momentum space at K{R) — R'^" (or higher) then the physically 
relevant predictions are not affected. 

The 0K41 theory, see [LL71], assigns to kq the value 3/4. Therefore the flows of physical 
interest should be described by (1.1) truncated at | A: | < K{R), i.e. 

ikk+iR X! Tik^-knkUk^_ = -k'^Uk + '£t^ (2.1) 

\kA<K{R) 
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where u{x)'^= X^/,-^o ^'~~^fc ^^'^ is the projection orthogonal to k; k = 2nn 
with n an integer components vector. 

Equation (2.1) admits an a prion bound on the energy E/2 = — X]^. ^^|wfeP + 
J2k^k ' y^-k which implies that asymptotically in time the energy is bounded by < 

2\W\fl{2ir)\ 

We shall call jiR the statistics^' of the NS equation defined as the probability distribu- 
tion on phase space {i.e. on the space of the velocity fields {w fe}, \ k\ < K{R)) which 
describes, at Reynolds number R, the stationary state averages of observables F{u), 
i.e. the probability distribution such that 

lim T-^ [ F{Stu)dt= [ F{v)nR{dv) (2.2) 

for almost all initial data u. The distribution jjji exists and is unique because of the 
chaotic hypothesis and it is also called the SRB distribution of the stationary state of 
(2.2). 

3. NS and GNS equations: viscosity and vorticity ensembles. Equivalence. 

For the purposes of a conceptual analysis stressing the analogy with the theory of 
ensembles in statistical mechanics we temporarily introduce a control parameter A in 
front of the Laplacian in (1.1) and in front of the — k'^ u k term in (2.1): bearing in 
mind, however, that we are interested in A = 1 only. The stationary distributions will 
then depend on A, R and will be denoted /Ua,_r so that the previously introduced SRB 
distribution /z_r, (1.1), is fiR = ^i^r with the new notations. 

The collection S^s of all the probability distributions hx^r is the collection of all the 
stationary states of the fluid, at varying Reynolds number R. It is an ensemble in the 
sense of statistical mechanics and it will be called the "viscosity ensemble" for the NS 
equations. 

The second idea is that the same fluid can be studied, rather than by the NS equations 

(2.1), by considering the Euler equations subject to a dissipation mechanism that keeps 

the vorticity 3"^= jyidu)"^ dx boimded. This "thcrmostatting" effect can be achieved 
by imposing various types of forces th on the system so that <S = ccmst leading to 

u+Ru-du=-dp+ ip + th, d-u = 0, R = fL^i/-^ (3.1) 

As an example we can consider the force obtained by imposing the constraint that <S 
remains identically constant via Gauss' principle of least effort, c.f.r. [Ga96b], appendix. 
It corresponds to, c.f.r. [Ga96b] 

u + Ru ■ d u = — dp+ (fi + vciu) Au, d-u = (3.2) 

where i^ailL) is an easily determined multiplier defined so that S is exactly conserved, 
namely 

Jy {ip -Am — RAu ■ {u ■ du)) dx 

^oiy) = ^ — fyiAuydx ~ 

where V is the container region. 

The collection £ of the statistics jls,R for (3.1) will be called the "vorticity ensemble" . 
We establish a correspondence between elements of the ensembles £ and £ by saying that 
two elements hx^r G £ and i^s,R & £ are correspondent if 

I^sM'^g) = A (3.4) 
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Wc shall call ^HocaP an observable F{u) that depends only on finitely many Fourier 
components of the velocity field u (this is locality in momentum space) and C is the 
family of the local observables. The analogy with statistical mechanics is quite manifest 
if the following conjecture holds, [Ga95b], 

Equivalence conjecture: Let I-1\.r G £ and ^s,R 'E £ be corresponding elem,ents of the 
viscosity and vorticity ensembles, i.e. if S and A are related by (3.4), then it is 

hm — = 1 (3.5) 
for all local observables F £ £ with non zero average. 

In other words the statistics of the irreversible NS equation (1.1) and of the reversible 
"Gaussian NS equation" (3.1), called GNS equation, form two equivalent ensembles of 
stationary states of the fluid. By "reversible" we mean that there is an involutory map 
/ of phase space which anticommutes with time evolution, i.e. 

1^ = 1, ISt = S-tI for all t (3.6) 

and the GNS equations are reversible because i^g{u) is odd under the transformation 
I u{x) = — u{x). A similar conjecture has been proposed in certain models of non 
equilibrium statistical mechanics, [Ga96b], [Ru99b], [Ga99d]. 

For reference purposes we write explicitly the GNS equations with the 0K41 cut-off 
ikk+iR XI Tik^-k^kUk^ = -yG{u)k'^Uk + <£^ (3.7) 

\'k.\<K(R) 

with the same notations of (2.1). 

The analogy with equilibrium theory of ensembles is: the parameter R plays the role 
of the volume while A that of the temperature and S that of the energy. Therefore the 

viscosity ensemble is the analogue of the canonical ensemble and the vorticity ensem- 
ble is the analogue of the microcanonical ensemble. The i? — »■ oo is analogous to the 
"thermodynamic limif , [Ga99a]. 

We see also why it is useful to introduce the parameter A: if we stick to A = 1 then 
effectively we consider only a single stationary state jJi^x = hr and not an ensemble: 
this state is "the same" (in the sense (3.5)) as the state fiSR,R ff Sr is so defined that 
jj'SR,R{'^G) = 1- The parameter A will be set to its physical value 1 from now on. 

The conjecture of equivalence was proposed in [Ga97a] and discussed in several other 
papers, see for instance [Ga97b]. It has been investigated by simulations in [RS99] with 
results that seem moderately satisfactory. 

§4. Time reversal and fluctuation theorem. 

We now consider the NS equation (3.7) and try to find some of its properties under 
the assumption that it is equivalent to the corresponding GNS equation, i.e. (3.7) with 
~k.'^ Ikk replaced by —VG{u)k'^Uk - assume the chaotic hypothesis and the 0K41 
cut-off and furthermore that 

Transitivity and axiom C: Either the full ellipsoid in phase space 

{u\ E k^\^k? = SR} (4.1) 

\k\<R'^o 
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is densely visited by the evolutions of all data starting on it a,part from a zero volume set, 
or alternatively the evolution on this ellipsoid verifies a geometric property called "axiom 
C". 

Axiom C says that if the system is not transitive because there is an attracting set 
A that is sm,aller than the fuU phase space {i.e. the elhpsoid (4.1) in this case) then, 
considering the simple case in which tliis happens because in phase space there are just 
a non dense attracting set and a repeUing set (also not dense), 

(1) the attracting and the repelhng sets are smooth manifolds and all their points, but a 
set of zero surface area, generate dense trajectories on them, and 

(2) the stable manifold of the points on the attracting set crosses transversally the re- 
pelling set and viceversa the unstable manifold of a point on the repelling set crosses 
transversally the attracting set. 

for details, which will not be really necessary here, we refer to [BG98]. This implies that 
either the system is transitive or that its restriction to the attracting set is transitive. 
The interest of the Axiom C notion is that it is a geometric property that has a re- 
markable consequence for systems admitting a time reversal symmetry I but with an 
attracting set A which is not the full phase space and, therefore, is mapped by I onto 
a repelling set I A different from A. If the axiom holds one can define, [BG98], a map 
P : A < — > / A, of the attracting set A on the repelling set I A which commutes with time 
evolution and with I and 



i.e. the restriction of the transitive evolution St to the attracting set A is still reversible, 
although it is such for a new tim,e reversal operation, namely I P, see [BG98] , [GaQSb] . 

// a reversible evolution verifies axiom C and depends on a parameter and, as the pa- 
rameter varies, it develops an attracting set A =/= I A that is not the full phase space then 
the restriction of the evolution to the attracting set is time reversible with respect to a 
new time reversal symmetry. 

In other words in axiom C systems time reversal symmetry / cannot be really broken: 
if there is a spontaneous breakdown (such has to be considered the "breaking", as a 
parameter varies, of phase space into an attracting set A smaller than phase space and 
a repelling set I A different from I A, [GaQSb]) a new time reversal PI is "spawned". 

The axiom C property is stable under perturbations: changing slightly parameters a 
system keeps this property if it has it to begin with. The transitivity (or axiom C) 
property is relevant because of the following theorem 

Theorem (fluctuation theorem): Let —a{u) be the divergence of the GNS equations 
(3.7) and let —{cr)+ be its stationary average with respect to the SRB distribution. 
Then the (dimensionless) quantity 



which we call "average over a time span t of the (cUmcusionlcss) phase space contrac- 
tion at u " has a probability of being in the interval \p,p + dp] of the form TTrip) dp = 
const e'^^P^^+^<^^^ with 



IPSt = S^tIP 



(4.2) 




(4.3) 



Q{-p) = ap)-{^)+p 



(4.4) 



for all p. 
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This theorem can be found in [GC95] for evolution maps and in [Ge98] for flows (which 
is the version we use here): see also [Ru99a]. The quantity p depends on u. The quantity 
(cr)+ is also called "average entropy creation rate" and p = p{u) is the dimensionless 
entropy creation rate averaged over a time r and in the point u : sec [An82] , [Ru96] , 
[GR97], [Ru99a]. We recall that entropy in systems out of equilibrium is not defined 
(yet) so that this name needs not be taken too seriously and it might eventually reveal 
itself inappropriate. 

The above result should not be confused (as it is conceptually and technically different) 
with other apparently similar statements, see [CG99]. It was discovered as an experimen- 
tal relation in a nmnerical simulation, [ECM93] , where the role of the SRB distributions 
and of time reversal were also suggested to be a possible reason for its validity: this led 
to its proof for Anosov maps in [GC95] and for flows in [Ge97] . 

A key feature of the theorem is that it contains no free parameters: its generality 
makes it a mechanical identity in the same sense, although of course of not comparable 
importance, as the heat theorem of Boltzmann, [B066] , [Bo84] , see also [Ga99a] . 

In the case of axiom C systems (4.4) still holds, because the evolution restricted to A 
is transitive and reversible by (4.2), hut a lias to be replaced by the contraction rate ctq 
of the surface area of the attracting set A. The quantities a and ctq seem unrelated; 
however there are important cases in which the total phase space contraction a and the 
contraction of the surface element of the attracting set (Tq a-re proportional: ctq = "dga with 
i?o a constant factor (or varying on a slower time scale than a itself). Then (fTo)+ = i?(f)+ 
with 1? = {'do)+ and (4.4) becomes 

C(-p)=C(p)-t^H+ (4.5) 

for all p. The GNS case is not among the (important) cases in which a and ctq are 
proportional, see [DM96], [WL98],[BGG97] and [BG97], although heuristic arguments 
can be given, [Ga97a], suggesting that nevertheless a relation like (4.5) might hold. 
In some cases in which proportionality between a and ctq can be established, at least 
on heuristic grounds, the proportionality factor is just 1 — d{A)/d if d is the dimension 
of phase space and d{A) is the dimension of the attractor, but unfortunately such cases 
are very rare, [BGG97], [BG97]. 

Finally it is worth writing explicitly the expression of the phase space contraction (t(m ) 
for the equation (3.7) 

a{u)^i k^)iyG{u)-{ [ Aip-Audx){ f [{Auf- 

\k\<K{R) "^^ 

-RAu- (A(u -du)) - R{Au) ■ {Au) ■ (du) - RAu ■ {Adu)u+ (4.6) 



+ i^{u)Au ■ A'^ujdx) I / {Auf dx 

Jv 



In this expression (straightforwardly derived by imposing that S is exactly constant on 
motions verifying (3.7)) the flrst term seems to be the dominant one at large R so that 

a(u)~( ^ k'^)y{u)'^= 1C{R)vg{u), IC{R) (x R^''°+^ (4.7) 

\k\<K{R) 

which, if the side Lq of the box is not L = 1 would be written with IC{R) replaced by 
K^LoiR) = T,\k\<K{R) with k = 27rn/Lo. 
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§5. Fluctuation patterns and an extension of Onsager— Machlup fluctuations 
theory. 

A physical interpretation of the flnctuation theorem, when it holds, can be fomid along 
with proposals for its test in experiments. We need first some consequences of the 
(technique of proof of the) fluctuation theorem. 

Given an observable H{u) we say that in its evolution observed over a time interval of 
size r it follows a pattern t h{t) if F{Stu) = h{t) for t e [— t/2, t/2]. We assume that 
F has well defined time reversal parity £ = ±1: F{I u) = £ F{u), for simplicity; and we 
say that the pattern Ih{t) = eh{—t) is the time reversed pattern of h. 

Fluctuation patterns are the main object of analysis in the theory of Onsager Machlup 
which deals with the probability of observing a fluctuation pattern h for an observable 
in the linear response regime {i.e. strictly speaking it deals with derivatives of various 
quantities with respect to the strength of the forcing terms evaluated at forcing). 

The following theorem can be regarded a result of the same type without the restriction 
that the system is in the linear response regime. 

Theorem (entropy creation as a fluctuations driver): 

Consider a time reversible evolution verifying the chaotic hypothesis and transitivity. 
Let H,K be two local observables (of given time reversal parity) and denote ij.r,t,p the 
SRB distribution conditioned to a (dimensionless) phase space average contraction p 
over a time span r. Let h, k be two fluctuation patterns for H, K and let Ik, Ik he their 
time reversal patterns. Then if HR^r,pipsitteTn oi H = h) denotes the probability that H 
follows the pattern h in the time span r in which the average dimensionless phase space 
contraction is p it is 

A* ij,r,p (pattern oi H = h) (pattern of H = Ih) 

Atfl,r,p (pattern of K = A;) /iK,,-,-p (pattern of K = Ik) 

for large t. If the system verifies axiom C the same result holds with IP (see (4-2) 
replacing I (without requiring any relation between the total phase space contraction rate 
and the rate of contraction of the surface of the attractor). 

In other words the relative probability of fluctuation patterns of H and K observed in 
a time span t during which the average entropy creation rate is p {a)+ are the same as 
those of the time reversed patterns in a time span r in which the average entropy creation 
rate is the opposite: —p (cr) + . 

This allows us to give a physical interpretation to p: namely if we look at the evolution 
on time laps of size r we see that the average entropy creation rate p will be usually 
p = I and the probability of observing p ^ 1 will be rare and the fraction of times we 
shall observe it is e*^'»*^^^~^^'^^^'^: hence events in which p ^ 1 will be rare and random 
{i.e. intermittent) and they will take place at rate ({!) — C{p)- 

The above theorem shows that when p is significantly different from 1 ^Hhings go very 
wroncf^ . The frequency of findings of a time interval of size r during which the time 
reversed patterns are relatively as probable as the normal patterns will be given by 
^-('^)+T no matter which observable H we look at: an independence property that can 

be checked, in principle, in an experiment. 

Hence a physical interpretation of p is that it is a quantitative measurement of the degree 
of reversibility that is observed. The larger 1 — p is the more ^^unintuitive behavior'^ 
will be observed. For p = —1 everything would be dramatically different from what 
expected.-^ The time intervals during which anomalous behavior is observed are rare so 



"If entropy creation rate could be changed in sign for a minute around Niagara falls then during that 
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that their manifestation is intermittent and we call this phenomenon entropy driven 
intermittenci/' : the function (^(p) describes the phenomenon quantitatively. 

§6. Entropy driven intermittency. Observability. 

We now address the question: "is this intermittency observable"? is its rate function 
C(p) measurable? 

Clearly a+ and ({p) will grow with the size of the system i.e. with the number of degrees 

of freedom, at least, which ^5 > 00 so that there should be serious doubts about the 

observability of so rare fluctuations. 

However if we look at a small subsystem in a little volume Vq of linear size Lq we can 
regard it again as a fluid enclosed in a box Vq described by the same reversible GNS 
equations. We can imagine, therefore, that this small system also verifies a fluctuation 
relation in the sense that if, c.f.r. (4.7), (3.3) 

J^^^ (£ ■ An - RAti ■ (11 ■ diiY) dx (g.l) 



i/g(w) 



then it should be that the fluctuations of a averaged over a time span r are controlled 
by rate functions (v {p) and Cvb (p) that we can expect to be, for R large 

Cv{p) = ap)lCL{R), and 

Cvoip) = ap)ICLo{R), (6-2) 

(crVo)+ = CT+/Cz,„(i?) 

We recall that, c.f.r. (4.7), JClq (R) ^J2\k \<k{Rlo) — ^ where Rl^ is the Reynolds number 
on scale Lq which from the 0K41 theory is Rlo — {Lo/L^^^R. So that K-LaiR) oc 
{Lq/L)^R^^^^. Hence if we consider observables dependent on what happens inside Vq 
and if Lq is small so that ICl„{R) is not too large and we observe them in time intervals 
of size T then the time frequency during which we can observe a deviation "of size" 1—p 
from irreversibility will be small of the order of 

e(C(p)-C(i))T/Cz,o(ii) (g_3) 
for r large, where the local fluctuation rate C(p) verifles (assuming transitivity or axiom 

a-p) = ap)-a+pi} (6.4) 

with 1? = 1 in the transitive case and perhaps ^ 1 when the attracting set is smaller than 
phase space. 

Therefore by observing the frequency of intermittency one can gain some access to the 
function ((p)- 

Note that one will necessarily observe a given fluctuation somewhere in the fluid if Lq 
is taken small enough: in fact the entropy driven intermittency takes place not only in 
time but also in space. Thus we shall observe inside a box of size Lq "somewhere" in the 
total volume V of the system a fluctuation of size 1—p with high probability if 

(i/Lo)'e«f)-^(i)) "'=^o(«) ~ 1 (6.5) 



minute their water would be more likely to go up rather than down". One "just" has to change the 
sign of the entropy creation rate, no extra eflfort needed! 
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and the special event p = —1 will occur with high probability if 

(L/Lo)^e-^+^'=^o(«) ~ 1 (6.6) 

by (6.4). Once this event is realized the fluctuation patterns will have relative probabil- 
ities as described in §5. 

An attempt at interpreting the experiment performed by Ciliberto and Laroche on 
convecting water at room temperature in terms of the above theory is in [Ga99d] . 

The idea and the possibility of local fluctuation theorems has been developed and tested 
first numerically, [GP99], and then theoretically, [Ga99c], by showing that it indeed works 
at least in some models (with homogeneous dissipation like the GNS and NS equations) 
which are simple enough to allow us to build a mathematically complete theory of the 
phase space contraction fluctuations. 

Of course if the quantity {<Jv)+ and p {(7v)+ could be measurable, like the equilibrium 
entropy, in terms of heat ceded to the various thermostats acting on the system divided 
by their temperature, then the whole theory could be even more easily subjected to 
experimental check as we could directly measure the rate function Qip) at least in small 
(but still macroscopically large) subvolumes. But the interpretation of the phase space 
contraction as a "physical entropy" (a concept that, however, we mentioned as still 
requiring a definition in non equilibrium physics) is quite controversial, e.g. see [Ho99] 
p. 236 and p. 240, as any statement about entropy is doomed to be. 

I adhere to the point of view, [An82], [GR97], [Ru99a], that the right definition of 
entropy creation rate in systems out of equilibrium is just the phase space contraction 
rate: but the connection with measurable entities of the quantity so defined is (therefore) 
an open problem. Nevertheless we have seen that there are already quite a few checks 
to test the theory that are already possible (and necessary given the large number of 
assumptions that must be made to obtain them). 

Acknowledgements: Work partially supported by Rutgers University and MPI via a 
grant # ?????. 
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